Proplyds and Massive Disks in the Orion Nebula Cluster Imaged 

with CARMA and SMA 



J. A. Eisner 1, , R. L. Plambeck 

University of California at Berkeley 
Department of Astronomy 
601 Campbell Hall 
Berkeley, CA 94720 

j ae@astro . berkeley . edu , plambeck@astro . berkeley . edu 

John M. Carpenter, S. A. Corder 

California Institute of Technology 
Department of Astronomy MC 105-24 
Pasadena, CA 91125 

jmc@astro . caltech. edu , sacOastro . caltech. edu 

and 

C. Qi, D. Wilner 

Harvard- Smithsonian Center for Astrophysics 
60 Garden Street, Mail Stop 42 
Cambridge, MA 02138 

cqi@cf a . harvard . edu , dwilner@cf a . harvard . edu 



ABSTRACT 



We imaged a 2' x 2' region of the Orion Nebula cluster in 1.3 mm wave- 
length continuum emission with the recently commissioned Combined Array for 
Research in Millimeter Astronomy (CARMA) and with the Submillimeter Array 



1 Miller Institute for Basic Research in Science, Berkeley, CA 94720 
2 Steward Observatory, University of Arizona, Tucson, AZ 85721 



-2- 



Our mosaics include > 250 known near-IR cluster members, of which 36 
are so-called "proplyds" that have been imaged previously with the Hubble Space 
Telescope. We detected 40 sources in 1 mm continuum emission (one of which is 
the BN Object), and several of them are spatially resolved with our observations. 
33 detected sources are known near-IR cluster members, of which 11 are proplyds. 
The 1 mm emission from the majority of detected sources appears to trace warm 
circumstellar dust. However, for many of the proplyds, which are located close to 
the Trapezium stars, the millimeter wavelength fluxes are dominated by thermal 
free- free emission from hot, ionized gas. Dust masses inferred for detected sources 
range from 0.01 to 0.5 M . For the ~ 225 known near-IR cluster members not de- 
tected in our 1 mm observations, images toward the positions of near-IR sources 
were stacked to constrain the mean 1 mm flux of the ensemble. The average flux 
is detected at the > 4a confidence level, and implies an average disk mass of 
~ 0.001 M Q , approximately an order of magnitude smaller than the minimum 
mass solar nebula. Most stars in the ONC thus do not appear to currently pos- 
sess sufficient mass in small dust grains to form Jupiter-mass (or larger) planets. 
Comparison with previous results for younger and older regions indicates that 
massive disks evolve significantly on ~Myr timescales. We also show that the 
percentage of stars in Orion surrounded by disks more massive than ~ 0.01 M 
is substantially lower than in Taurus, indicating that environment has an impact 
on the disk mass distribution. Disks in Orion may be truncated through photo- 
evaporation caused by the intense radiation field of the Trapezium stars, and we 
see marginal evidence for such a scenario in the spatial distribution of massive 
disks within the cluster. Our data show no statistically significant correlation 
between disk and stellar masses, although we see hints of a higher percentage of 
massive disks around lower mass stars. 

Subject headings: Galaxy: Open Clusters and Associationsdndividual: Orion, 
Stars: Planetary Systems:Protoplanetary Disks, Stars: Pre-Main-Sequence 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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Introduction 



The existence of protoplanetary disks around young stars is now firmly established. 
High resolution images from optical to radio wavelengths have shown disk-like morphologies 
and Keplerian rotation profile s around a number of young stars (e.g.. lO'Dell &: WoneJll99d 



1995; 



McCaughrean fc O'Delll 1999: IPadeett et al l 11999 : Eisner et all 12004 ; iKoerner fc Sargent 



Dutrev et all Il996l : Iwilner et al.ll200oh . Moreover, observations of near-IR emission 



from young stars in excess of that expected from their stellar photospheres imp ly that most 



stars aged less than a fe w million years possess inner circumstellar disks (e.g.. IStrom et al. 
19891 : lHaisch et al.ll200lh . 



Protoplanetary disks are the birth-sites of planetary systems, and the ubiquity, proper- 
ties, and lifetimes of disks constrain the timescales and mechanisms of planet formation. The 
mass distribution of protoplanetary disks is especially important since disk mass is related 
to the mass of planets that may potentially form. For our own solar system, the masses of 
planets and other bodies can be used to reconstruct a minimum-mass solar nebula (MMSN) 
describing the amount of solar-composition material needed to build the solar system. De- 
pending primarily on the core masses (and hence chemical compositions) of Jupiter and 
Saturn (which are not known precisely), estimates of the MMSN range from ~ 0.01-0.1 M 
(IWeidenschiUingi 119771 ). Such disk masses are also required by planet forma tion models to 



build giant planets on timescales shorter than inferred disk lifetimes (e.g., lHayashil 11981 



Alibert et al.ll2005l ). The MMSN is thus an informative benchmark against which to gauge 
the potential of disks around other stars to form solar systems like our own. 



A widely used method (e.g., iBeckwith et al.lll990l ) for measuring disk masses is to ob- 



serve emission from optically thin dust, and then use assumed dust grain properties to 
convert observed fluxes into dust masses. An assumed gas-to-dust ratio is then used to 
estimate the total (gas+dust) circumstellar mass. At short wavelengths (A < 10 /mi), the 
dust in protoplanetary disks is optically thick even for masses < 10~ 6 M . Observations 
at sub-mm and mm wavelengths are necessary to measure optically thin dust emission, and 
hence to determine the total mass of dust in the disk. 

Several investigators have carried out comprehensive single-dish mm and sub-mm con- 



(Bee 


kwith et al.lll990: 


Osterloh & Beckwith 1995 


: Motte & Andre S 


J001; 


Andrews & Williams 


2005 


). o ODhiuchi (Andre & Montmerle! Il994l: I 


^uernb 


erger et al. 


Il99< 


3; 


Motte et all Il998l: 


Andrews & Williams! 


20071). Lupus (Nuernberger et al. 


1997). Chamaeleon I (IHenning et al. 


1993 


). Seroens (Testi k Sargent 1998). and MBM 12 


(Itoh et al 


2003; 


Hogerheiide et al. 


2002 


). In Taurus, 22% of stars surveyed appear to possess disks more massive than ~ 0.01 



M , and the median disk mass is 5 x 10~ 3 M ([Andrews fc Williamsl 120051 ) . The fraction of 
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mass i ye diskaj and the median d isk mass is comparable in p Ophiuchi ( Andre fc Montmerle 
19941 ; [Andrews fc William] l2007h . 



Low-density star forming regions are not the typical birth-site s of stars; rather, most 



stars form in rich clusters like the Orion Nebula (ILada et al.l Il99ll . 1 19931 ; ICarpenterl 12000 



Lada fc Ladall2003l ). Isotopic abundances in our s olar system suggest that it, too, may have 



form ed in a dense, Orion-like environment (e.g., iHester fc Deschl 120051 : IWilliams fc Gaidos 



20071 ). Expanding millimeter continuum surveys to include rich clusters allows the determi- 
nation of the frequency and evolution of massive disks in typical star (and planet) formation 
environments. The high stellar density in rich clusters also allows the assembly of good 
statistics since many disks can be mapped at once. 

The main challenge to observing rich clusters at (sub-)mm wavelengths is that very 
high angular resolution is required to resolve individual sources and to distinguish compact 
disk emission from the more extended emission of the molecular clouds in which young 
clusters are typically embedded. Single-aperture mm-wavelength telescopes lack sufficient 
angular resolution, and to date, only three rich cl usters have been observed with mm 



wavelength interferometers: the Orion Nebula cl uster (IMundv et al 



Williams et al 



2005 



( lEisner &: Carpenter! 120031 ). 



Eisner fc Carpenter! 120061 ) . IC 348 (jCarpenter 



9951: bally et alill998bl : 
2002( 1. and NGC 2024 



These interfere metric surveys of rich clusters have detected very few disks with > 0.01- 
0.1 M of material, in large part because of limited sensitivity and areal coverage. The 
most recent observations of Orion detected emission from several massive (> 0.01 M ) disks 



( Williams et al.ll2005l ; lEisner fc Carpenter |2006), while upper lim its from other surveys range 



from ~ 0.025-0.17 M (IMundv et al.lll995l ; iBally et al.l Il998bl ). Considering as ensembles 
the large numbers (> 100) of young stars included in the cluster surveys allowed estimates 
of mean disk masses of ~ 0.002, 0.005, and 0.005 M^in IC 348, NGC 2024, and the ONC, 



respectively (ICarpenterl |2000| ; lEisner &: Carpenter! 120031 . 120061 ). Thus, it appears that many 
stars aged < 1 Myr still possess massive circumstellar disks, although more sensitive obser- 
vations are necessary to detect directly large numbers of massive disks at a range of ages, 
and thereby constrain the mass distribution and evolutionary timescales. 

Here we present a new 1.3 mm wavelength interferometric survey of the Orion Nebula 
cluster (ONC), a young, embedded stellar cluster that includes the bright, massive Trapez- 
ium stars. Our observations make use of the Submillimeter Array (SMA) and the recently 
commissioned Combined Array for Research in Millimeter Astronomy (CARMA). The com- 



2 Here and throughout the text, "massive disks" refer to disks with mass comparable to or greater than 
0.01 Mq, the lower range of estimates for the minimum-mass solar nebula. 
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bination of data from these two instruments yields a map of the mm-wavelength continuum 
emission in Orion with unprecedented sensitivity, angular resolution, and image fidelity. 

The Trapezium region contains hundred s of stars within a several arcm inute radius, and 
pre-main-sequence evolutionary models (e.g.. iD'Antona fc Mazzitellilll994j ) fitted to spectro- 
scopic and/or photo metric data indicate that most stars a re less than approximately one 



million years old (e.g., Prosser et al. 1994 ; Hillenbrand Il997 ). The standard deviation in the 



distribution of inferred stellar ages is < 1 Myr ([Hillenbrand! 1 199 71 ) . Our observations thus 
provide a snapshot of millimeter emission around a large number of roughly coeval young 
stars. 

With the large number of stars in the ONC, we can also investigate the correlation of disk 
properties with stellar and/or environmental properties. Previous investigations of near-IR 
excess emission showed the inner disk fraction for stars in Orion to be largely independent 
of stellar age and mass, although there are indication s of a paucity of disks around very 
massive stars ([Hillenbrand et al.lll998; lLada et al.ll2000n. In a ddition, the inner disk fraction 
may decrease at larger cluster radii ([Hillenbrand et al.lll998l ). The millimeter observations 
presented here enable investigation of how such stellar and environmental properties correlate 
with disk mass. 

We adopt a distance to the ONC of 400 pc, based on recent trig onometric parallax 
measurements of several st ars (jSandstrom et al.ll2007l ; iMenten et al.ll2007l ) and orbital fitting 



for a spectroscopic binary (IKraus et al 



20071 ). This is substantially lo wer than the yalue o f 



480 pc computed based on statistical parallax of water maser spots (IGenzel et al.l 198ll ). 
which was adopted in previous studies of the ONC (e.g., Eisner fc Carpenterl 120061 ). We 
discuss below the importance of this revised distance on our results. 



2. Observations and Data Reduction 

2.1. CARMA Observations and Calibration 

We mosaicked a 2' x 2' region toward the ONC in A1.3 mm continuum with CARMA 
between October and December, 2007. CARMA consists of six 10-m antennas and nine 6-m 
antennas situated at 2200 meters elevation at Cedar Flat in the Inyo Mountains of California. 
With a total of 15 antennas, CARMA provides 105 baselines, enabling excellent coverage of 
the uv plane and hence high image fidelity. Two different array configurations ('C and 'B') 
were used to obtain antenna separations ranging from 30 to 946 meters. 

Continuum data were recorded in six ~ 500 MHz bands covering the frequency ranges 



- 6- 



221.75-223.25 GHz and 227.25-228.75 GHz from the receivers' lower and upper sidebands, 
respectively. Each band consists of 15 channels. Spectral line emission in Orion is mostly 
resolved out by these observations; that is, across most of the mapped region it is detected 
only weakly because it is spatially extended relative to the interferometer fringe spacings. 
The spectral lines most visible on the shortest baselines, mostly toward the BN/KL region 
in the northwest corner of the mosaic, are the 11(1,11)-10(0,10) transition of SO2 at 221.965 
GHz and the 25-24 transition of HC3N at 227.402 GHz. Even these lines are almost com- 
pletely resolved out for projected antenna spacings > 70 kA (used to generate the final maps; 
see §2.3p . We therefore assume our bands are effectively line-free. 

The mosaic consists of 16 pointing centers (Figured]), separated by 26". This separation 
is comparable to the FWHM beamwidth of the 10-m antennas, but \/3 smaller than the 
FWHM beamwidth of the 6-m antennas. A 2-D mosaic is Nyquist sampled if pointings are 
separated by <FWHM/V3. Thus the CARMA mosaic is Nyquist sampled only for the 6-m 
dishes. However, Nyquist sampling is not crucial since we are interested only in compact 
sources, rather than extended emission. Simulated CARMA mosaics of a synthetic star 
field showed that 26" spacings provided the best balance of sensitivity and areal coverage to 
maximize the number of sources detected. 

Each night we interleaved 16-minute observations of the ONC mosaic, with 1 minute 
integration time on each pointing center, with 3-minute observations of the phase and am- 
plitude calibrator, J0530+135. Observing the mosaic in its entirety every 20 minutes and 
repeating this multiple times during the night ensures a high quality synthesized beam and 
equal sensitivity for each pointing center. The total integration time for the maps was 34 
minutes per pointing center. 

Telescope pointing was checked and updated on 20 minute intervals using optical coun- 
terparts very near to the source. These objects were observed with optical cameras mounted 
on the dishes, and radio-optical offsets were calibrated periodically. Gain stability, especially 
in 1 mm observations like those presented here, is found to be enhanced through this method 
(Corder, Carpenter & Wright, in prep). 

We measured a 1.3 mm flux density of 3.0 ± 0.3 Jy for J0530+135 using Uranus as a 
primary flux calibrator, based on observations during several nights in October when both 
sources were observed. J0530+135 was also used to calibr ate the passband. All calibrations 



for these data were performed with the MIRIAD package (ISault et al.lll995l ). 
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2.2. SMA Observations and Calibration 

We mosaicked a triangular region approximately V.7 on a side, consisting of three point- 
ings, with the SMA between September, 2005, and February, 2006 (Figured]). The SMA 
consists of eight 6-m dishes near the summit of Mauna Kea in Hawaii. However, for our 
observations only seven antennas were available, providing 21 baselines between 10 and 220 
m. The uv coverage for our SMA observations is substantially sparser than for our CARMA 
observations, and thus the image fidelity is worse (which means, for example, that strong 
emission is scattered more strongly into other regions of the map). However, the high alti- 
tude of the SMA enables very low opacity observing conditions. The low opacities, combined 
with increased observing time per pointing, lead to substantially better sensitivity relative 
to the CARMA observations. 

Double sideband receivers were tuned to a local oscillator (LO) frequency of 225.333 
GHz. The SMA digital correlator is configured with 24 partially overlapping bands of 104 
MHz width in each sideband. Each sideband provides 2 GHz of bandwidth, centered ±5 
GHz away from the LO frequency. The double-sideband (DSB) system temperatures were 
between 80 and 200 K. 

As for the CARMA mosaic, we observed the mosaic multiple times throughout the night, 
obtaining equal integration times for each pointing position. The pointings in the mosaic 
were separated by ~ 44", the approximate FWHM for the SMA dishes at this observing 
wavelength. As for the CARMA map, the SMA mosaic is larger-than-Nyquist sampled, 
which provides enhanced areal coverage compared to a Nyquist-sampled mosaic. 

We used J0423-013 and 3C120 as gain calibrators for these tracks, with flux densities of 
1.4 and 1.0 Jy, respectively, derived using Uranus as a primary flux calibrator. We estimate 
that the absolute flux scale is uncertain by ~ 10%. We calibrated the passband using the 
quasar 3C454.3, and Uranus where available. Because the CO(2-l) transitions are present in 
the observing window, we edited out the parts of the band with strong lines and generated a 
line-free continuum channel. All calibrations were performed using the SMA adaptation of 
the IDL-based data reduction package MIR developed at Caltech; calibrated data were then 
converted into MIRIAD format for further processing. 

2.3. Mapping 

We made mosaics of our CARMA and SMA datasets individually, and after combining 
the two datasets in the uv plane. For the individual and combined datasets, we mosaicked 
the individual pointings into a single image, weighting the data by system temperature and 
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by uv distance (with a "robust" parameter of 0.5), then de-convolved and CLEANed (all 
using MIRIAD). The angular resolution afforded by the longest baseline data in our maps 
(from the CARMA B-array) is ~ 0''3. Our mosaics have 0''1 pixels, which ensures adequate 
sampling of individual resolution elements. 

Since we are primarily interested in compact disk emission, we eliminated uv spacings 
shorter than 70 kA (i.e., projected baselines shorter than 93 m) in order to reduce contam- 
ination from bright extended emission. The eliminated spacings correspond to size scales 
larger than ~ 3". The cutoff value was chosen to minimize the RMS background noise in 
the CLEANed images; we measured the RMS for uv cutoff radii of 50, 60, 70, and 80 kA, 
and found the 70 kA cutoff to be optimal. 

Mosaics produced from our robust-weighted data with r uv > 70 kA are shown in Figures 
[2HH In the Figures we have divided by the theoretical sensitivity at each location in the 
image, in order to visually down-weight the noisier edges of the mosaic (where there are fewer 
overlapping pointings); we do not divide by the sensitivity in the analysis presented below. 
We note that even within the uniform (theoretical) sensitivity region of the mosaic, the RMS 
varies substantially because of emission scattered from bright compact and extended sources 
in the BN/KL and OMC1-S regions located in the upper and lower right quadrants of the 
maps. 

For the SMA mosaic, the unit gain region (within which the theoretical sensitivity 
does not vary substantially) encompasses a roughly triangular region covering ~ 2 square 
arcminutes. The RMS of pixels within "clean" regions of the unit gain contour (i.e., away 
from the crowded BN/KL and OMC1-S regions) is 0.8 mJy The unit gain region of the 
CARMA mosaic encompasses a 2' x 2' area, with an RMS noise level (again, in clean regions 
of the map) of 2.3 mJy. For the combined map, which will be used for the bulk of our 
analysis, the unit gain region is slightly larger than for the CARMA-only mosaic, and the 
RMS noise level is 1.8 mJy. The synthesized beam has dimensions (FWHM) of 0''69 x 0"60 
at a position angle of 72°. 

3. Analysis and Results 

3.1. Detection Thresholds 

Because the map contains a large number of pixels, we must employ a fairly high detec- 
tion threshold to avoid random noise spikes if we search the image blindly. The mosaic area 
is approximately 35, 000 synthesized beams. With this number of independent pixels, one 
expects > 1 noise spike above the 4a level (assuming Gaussian noise). We therefore use a 
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5a detection limit, at which level <C 1 pixels are expected to show noise spikes. Because the 
noise varies greatly across the map, we calculate a locally in small sub-regions of the image. 

Specifically, a "local" a is computed in 10" x 10" (100 x 100 pixels) boxes around 
each pixel in the mosaic. For our detection thresholds to be meaningful, the noise must be 
well-characterized. However, poorly sampled extended emission leads to excess noise in the 
BN/KL and OMC1-S regions. Moreover, the noise increases toward the edges of the unit 
gain region because there are fewer overlapping mosaic pointings there. Detections in these 
areas should be treated with some caution. 

As a test of our detection threshold, we searched the maps for false detections below 
the —5a level. None were detected, confirming that 5a is a reasonable detection limit. In 
contrast, 28 sources were seen below the —4a level (most of them toward the edges of the 
mosaic or in the BN/KL and OMC1-S regions), demonstrating that 4a is not a sufficiently 
stringent detection threshold (and that the noise across our mosaic is not always Gaussian). 

Instead of blindly searching for detections, we can also use our prior knowledge of the 
locations of near-IR cluster members and search only these positions. For these ~ 250 
pre-determined positions, ~ 0.3 sources are expected to show emission above the 3a level 
from Gaussian noise. We can therefore try a 3a detection threshold, where a is the noise 
determined locally (as above) in 10" x 10" sub-regions centered on individual cluster member 
positions. Although the noise in the mosaic is not always Gaussian, this 3a threshold appears 
reasonable: none of the near-IR source positions were detected below the —3a level. 

Sources with 1 mm continuum emission at the > 3a level in our maps are deemed to 
coincide with near-IR cluster members if the mm peaks and near-IR source positions lie 
within 0"4 of each other. The estimated relative positional accuracy of 0"4 is the quadrature 
sum of uncertainties from centroiding the mm images (~ 0.5#bcam/sig na l-to-noise ~ 0"1), 
uncertainties in the absolute astrometry due to baseline errors (~ 0"2), and uncertainties in 
the near-IR source positions (~ 0"3). 

We detected 19 sources within the unit gain contour of our mosaic above the 5a level (Ta- 



bled]) . 12 of these are coincident with near-IR cluster members listed in lHillenbrand &: Carpenter 



( 120001 ). An additional 21 objects were detected above the 3a level toward positions of near-IR 
sources. 1 mm continuum images of sources detected in our mosaic are displayed in Figure 

El 



While the BN objec t is detected, we defer discussion of this high- mass, embedded object 
(e.g., iGezari et al.lll998l ; iPlambeck et al.l 119951 ) to a later paper that examines the BN/KL 
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region m detaii. In the remaining discussion, we focus our attention on the sources detected 
at both infrared and millimeter wavelengths. While most of these objects do not have known 
stellar masses, they are likely to be low-mass stars based on the stellar mass distribution 
computed statistically for the ONC as a whole ([Hillenbrand &: Carpenter! 120001 ) . 



3.2. Angular Sizes of Detected Objects 

For each source detected in our 1 mm mosaic, we fitted a 2-D elliptical Gaussian to the 
emission. The synthesized, clean, beam for the combined mosaic is a 2-D Gaussian with 
FWHM of 0"69 x 0'.'60 at a position angle (north of west) of 72°. At the assumed 400 pc 
distance to the ONC, the core of the synthesized beam has dimensions of 240 AU by 280 
AU (again at a PA of 72°). For simplicity, we approximate this as 240 AU in the East- West 
direction and 280 AU in the North-South direction. 

For sources detected at a signal-to-noise ratio of ~ 5, the statistical uncertainty in the 
fitted FWHM is ~ 10%. Baseline errors or phase noise in our mosaics can broaden the 
apparent source sizes, however. We assume that a source is resolved only if the major or 
minor axis of the fitted FWHM is 25% larger than that of the synthesized beam. Objects 
for which the fitted Gaussian FWHM is smaller are considered to be unresolved. 

Approximately 25% (9/39) of detected sources are spatially resolved in our images (Ta- 
ble El). An additional nine objects (all proplyds) have been spatially resolved with HST 



Jvicente fc Alvesl l2005h . Thus, angular sizes are available for ~ 50% of our sources. The 



inferred radii for resolved sources range from ~ 90 to ~ 220 AU. For unresolved objects we 
can say only that the emission is confined to radii smaller than9 ~ 100 AU . For a sample of 



134 p roplyds with sizes measured with HST, the mean disk radius is 71 AU (IVicente fc Alves 



20051 ). Since the mean disk diameter is ~ 1/2 the size of the linear resolution of our ob- 
servations, it is not surprising that most of the sources detected in the 1 mm mosaic are 
unresolved. 



3 Although Source I is detected as a strong, individual object in our CARMA B-array data, we do not 
detect it in our combined SMA+CARMA mosaic because of confusion with the hot core. Source I will also 
be discussed in the later paper. 

4 Larger sources, when convolved with the synthesized beam, would produce measured sizes > 25% broader 
than the beam. 
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3.3. Distinguishing Dust and Free-Free Emission 



Since we are interested in using our observations to constrain the mass of circumstellar 
dust around our sources, we must account for potential contributions to the observed fluxes 
from sources other than dust emission. Free-free emission arises in hot ionized gas, and 
in the ONC such conditions may exist either in HII regions around high- mass stars (e.g., 
Garav et al.1 Il987t IPlambeck et alJll995h or in the outer regions of disks or envelopes tha t 
are irradiated by the hot Trapezium stars (e.g.. lO'Dell et al.lll993l ; iHenney fc 0'Del]||l9 99h 
While some sources in th e ONC have shown flares of cyclotron emission (e.g.. iBower et al. 
20031 ; iFuruya et al.l 120031 ) . we expect that such rare events will not contribute significantly 
to the 1.3 mm fluxes, and we do not consider them here. 

Because the spectral shape of free-free radiation differs from that of thermal dust emis- 



sion, comparing 1 mm measurements with longer-wavelength data enab 
these components. We use long-wavel e ngth fluxes from the liter a ture ( 



Mundy et al.l Il995l ; IZapata et al.l 120041 ; Eisner fc Carpenter 



es us to distingui sh 



Felli et a 



200d: iForbrich et al.ll2008h . In 



1993 



as 



2005) for the few objects where 



addition, we use 880 /im fluxes measured by lWilliams et al.l ( 
these are available. 

For a freely expanding, fully ionized wind with constant M, such as we expect for 
proplyds, free-free emission will have the following spectrum^: 



F„s — 



v,tum 



,-0.1 
\0.6 



if v > 

5: ^turn 



if v < 



(i; 



Here, Vtum is the frequency above which the wind is optically thin at all radii. We include 
a derivation of th i s resu lt in the appendix, and an alternative derivation can be found in 
Wright fc Barlow! Jl975h . 



For v < u tUTn , the inner parts of the wind are optically thick to free-free radiation. If 
we adopt a simple model fo r proplyd winds where M = 10 -7 M yr" 1 with spherical wind 



velocity of 20 km s 1 (e.g., IHenney &: Q'Delll 119991 ). and T e = 10 4 K, we can estimate the 
size of the optically thick region. Using Equation IA51 we obtain x r ~i ~ 1 AU at Al mm and 



~ 30 AU at A10 cm. Even at 10 cm, this is smaller than the likely wind-launching regions 
for proplyds, and free-free emission from most proplyds is likely t o be fully optically thin . 



For the highest measured mass loss rates of ~ 10 6 M yr 1 (e.g.. IHenney fc O'Delll 119991 ). 



In a previous paper ( Eisner fc Carpente rll2006h . we assumed that free-free emission originated from static 
HII regions rather than from winds, as in the present work. This choice affects only the long-wavelength 
behavior of the free-free spectrum, and is relatively unimportant to our analysis. 
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the optically thick regions of proplyd winds are ~ 130 AU for A ~ 10 cm. We therefore do 
not expect to see a spectral turnover (Equation [1]) for wavelengths < 5 cm. 

Emission from cool dust is added to this free-free emission to obtain a model of the 
observed flux. We assume that 

i^dust = i 71 230GHz,dust(^/230 GHzf + « = F 230 GHz,dust (^/230 GHz) 3 , (2) 



for j3 = 1 (e.g., iBeckwith et al.lll990l ). Other values of (3 can not be ruled out based on 
our data for most objects, and (3 = 0, corresponding to emission from optically thick or 
large-grained dust, is typically compatible with the data. 

We estimate the relative contributions of dust and free-free emission by fitting this 
model, F v = F Vj q + F u ^ust, to our measured 1.3 mm fluxes and to 880 /mi, 3 mm, 3.6 mm, 
1.3 cm, 2 cm, 6 cm, and 20 cm fluxes from the literature. For comparison, we also fit a 
dust-only model, described by Equation [2J with (3 = 1 and (3 = 0. For objects detected 
at centimeter wavelengths, we fit the dust-l- free-free model to the > 4 flux measurements 
for each source, and thus we are able to determine the three free parameters of the model, 
^tum, -Ft/,turn, and -F^dust- For sources with < 3 flux measurements (i.e., those undetected in 
centimeter wavelength surveys), we fit only the dust emission model to the data. 

Given the noise level of previous centime ter observations covering the entire region of 
our 1 mm mosaic (< 0.3 mJv; belli et "aHll993al lbl) and the measured 1 mm fluxes for detected 
objects (> 10 mJy), sources undetected at cm wavelengths are probably dominated by dust 
emission. For a source with 10 mJy flux at 1 mm, a non-detection at 10 cm implies that 
< 0.2 mJy, or < 2% of the measured 1 mm flux is due to free- free emission. For simplicity, 
we attribute 100% of the 1 mm fluxes to dust emission for these objects. 

Fluxes, from sub-millimeter to radio wavelengths, and models are plotted in Figure [6j 
and the fluxes due to thermal dust emission are listed in Table [2j Uncertainties for these 
dust fluxes are given by the la uncertainties of the model fits. The majority of detected 
sources appear to be dominated by dust emission. However, for the subset of the sample 
seen in optical emission or absorption with HST (the proplyds at the top of Tabled]), the 1 
mm fluxes are dominated by free-free emission. This probably reflects the relative proximity 
of proplyds to the luminous Trapezium stars. 
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3.4. Estimating Circumstellar Dust Masses 



The mass of circumstellar dust is related to the component of the 1 mm continu u m flu x 
due to dust emission. Assuming the dust is optically thin, and following iHildebrandl (119831 ). 



M, 



dust 



ft v,dust B u ( Td us t ) 



(3) 



Here, v is the observed frequency, SV.dust is the observed flux due to cool dust, d is the distance 
to the source, ft^dust = fto^/^o)^ is the dust mass opacity, Td us t is the dust temperature, 
and B u i s the Planck func t ion. We assume d ~ 400 pc, k = 0.0002 cm 2 g _1 at 1.3 mm, 
(3 = 1.0 (IHildebrandl Il983l ; iBeckwith et al.l Il990h . and T dus t = 20 K (based on the average 



dust temperature inferred for Taurus; Andrews & Williams 2005; see also the discussion 
in Carpenter 2002; Williams et al. 2005). The dust mass can be converted into a total 
circumstellar mass by assuming the canonical gas-to-dust mass ratio: M C i rcum steiiar = ^dust x 
100. Column 3 of Table [2] lists the estimated circumstellar masses for detected objects. 

Uncertainties in the assumed values of these parameters (notably k) lead to uncertaintie s 
in the derived masses (in an absolute sense) of at least a factor of 3 (e.g. JPollack et al.lll994j ). 
which are not included in the uncertainties listed in the table. Values of ft^dust and Td us t 
may also vary across o ur sample. Since the cluster population in Orion is roughly co-eval 



[e.g.. iHillenbrandl 119971 ). such effects should be minimal. However, there is some spread in 
stellar masses, which may lead to some range in these parameters. For example, since some 
of the objects in Tabled] may be massive stars, the millimeter flux may contain contributions 
from dust hotter than the assumed 20 K; if T dus t = 30 K, then the computed dust masses 
would be lower by a factor of 1.6. For the predominantly low-mass cluster population 
( [Hillenbrand fc Carpenter! |2000| ). from which the sources listed in Table [1] are drawn, the 
assumed values for fi^dust and Td us t presumably do not vary much, and the masses predicted 
by Equation [3] are probably valid in a relative sense to within a factor of two. 



3.5. Constraints on Dust Optical Depth 



We perform a simple test to determine whether the optical depth (r) might become 
comparable to or larger than unity in the systems under study. Using disk sizes measured 
from either our data or HST data (IVicente fc Alvesll2005l ). or limits on disk sizes from our 
observations, we compute the emission expected from optically thick dust with a temperature 
of 20 K: 2 

-^disk \ • i i r> t i -^disk 



Fv,r>i — B u {T AvLSt )Ti 



d 



cosi ~ 112 mJy 



100 AU 



(4) 
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Here i?disk is the disk radius, i is the inclination and d is the distance. For simplicity, we 
take i — 0, which leads to an upper limit on the flux for an optically thick disk of radius 
-Rdisk- Fluxes expected for optically thick dust for our sample are listed in Table [2J 

For most objects, the fluxes (or upper limits) expected for optically thick dust are 
substantially higher than our measured fluxes (Table [2]). These disks must be either highly 
inclined or optically thin. Since it is unlikely that all of the disks detected in our observations 
are edge-on (especially since these would be the dimmest portion of a sample of randomly 
inclined, optically thick disks), we take this as evidence for optically thin material. Although 
many of the disk sizes, and hence the expected optic ally thick fluxes, are u pper limits, the 



mean radius for proplyds in Orion of ~ 70 AU (e.g., IVicente &: Alvesll2005l ) would produce 



an optically thick flux of ~ 50 mJy, still higher than the majority of our measured fluxes. 

We therefore believe that most disks detected in our observation are composed largely 
of optically thin dust, and the circumstellar dust masses computed in §3.41 are reasonable 
for these sources. There are a few among the sources detected only at > mm wavelengths 
(MM8, MM13, and MM21) for which the measured fluxes are comparable to or larger than 
the expected optically thick fluxes. The dust in these objects is either optically thick or 
hotter than 20 K, as may occur around higher-mass (proto-) stars. 



3.6. Stacking Analysis 

With the large number of young stars contained within our mosaic, we can enhance the 
effective sensitivity by considering the ensemble of ~ 225 sources not detected individually. 
For each known cluster member within the mosaic that is not detected above the 3a noise 
level, we make a 10" x 10" sub-image centered on the stellar position. We weight the sub- 
image by the local RMS (determined as described in §3.11) . sum all of the weighted images, 
and divide by the sum of the weigh ts to normalize. W e exclude any cluster members known 
to have radio-wavelength emission (belli et al.lll993al fbh. 



The weighted image is shown in Figure [7J The average flux for the ensemble of non- 
detected sources is 0.9 ± 0.2 mJy, with a significance of > 4a. The peak flux is centered on 
the mean position of the near-IR sources (within the positional uncertainties of ~ 0"4), and 
resembles the synthesized beam core, indicating that the average source is compact. 

Since the positional uncertainties are comparable to the half-width half-maximum of 
the synthesized beam, the average flux seen in Figure [7J may be slightly reduced because 
different sources in the ensemble do not lie exactly atop one another. Assuming the positional 
uncertainties are random and Gaussian-distributed, one would expect a reduction in the 
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measured average flux of ~ 35%. Correcting for the potential flux-smearing, one would 
obtain an average flux for the ensemble of 1.2 ± 0.2 mJy. We verify this by integrating the 
central region of the average image over a region with four times the area of the synthesized 
beam; as expected, we find an integrated flux of ~ 1.2 mJy. 

Low-level free-free emission may contaminate the average image, and hence bias the 
average flu x inferred for the ensemble. The la sensitivity in cm-wavelength surveys is ~ 0.3 



mJy (e.g., iFelli et al.lll993bl ). We argued in §3.31 that gaseous winds in the ONC are likely 
to be optically thin to free-free emission, and hence that F u ^ oc u~ 0,1 . Thus we would 
expect free- free emission to be no stronger than ~ 0.2 mJy at 1 mm wavelengths for sources 
undetected in cm- wavelength surveys. This is comparable to the la sensitivity in our 1.3 
mm average image. Since the average image is detected above the 4<r level, > 75% of the 
average flux comes from dust. 



Discussion 



4.1. Nature of Detected Sources 



The sub-arcsecond resolution of the CARMA observations is enough to marginally re- 
solve some of the detected sources in the ONC, and in principle we could observe flattened, 
disk-like geometries. For example, Figure [8] shows that 177-341 has a disk-like morphology 
aligned with the silhouette disk seen by HST. While only a few sources can be well-resolved 
with our observations, HST observations show that many of the observed proplyds appear 
disk-li ke (IMcCaughrean fc O'Delll Il996l ; iBally et al.l Il998al ) , some even exhibiting silhouette 
disks ( jBally et al.ll2000l ). For the proplyds and well-resolved mm sources, the 1 mm emission 
evidently arises from disk-like distributions. 

Mid-IR emission is also observed toward many of the sources detected at 1 mm. 82 



of so urces (all except the "MM" sour ces and HC 495) are also seen at 3.6 /im (ILada et al. 



2004 ). and 48% are seen at 11.7 /im (jSmith et al.ll2005l ). While 3.6 /im emission may trace 
stellar photospheres and/or infrared excess, the 11.7 /zm emission provides direct evidence 
for circumstellar material at least out to radii of a few AU. Thus, many detected sources 
(the majority, if the "MM" objects are excluded) have evidence for inner circumstellar disks. 

More generally, where 1 mm emission is detected toward known near-IR cluster mem- 
bers, the fact that the near-IR light is visible despite the high extinctions (Ay > 300) that 
one would derive based on the amount of material needed to produce the 1 mm emission (for 
spherically dist ributed material) implies that the dust lies in fla ttened, disk-like distributions 
(see also, e.g., 



Beckwith et al.lll990l : Eisner fc Carpenter! 120031 ). 
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It is interesting to speculate as to the nature of sources detected at HKL bands and 
at mm wavelengths, but not at 11.7 /im. 13/32 (40%) of mm and near-IR detected sources 
fall into this category. It is possible that some of these are transitional disks. The HKL 
emission may trace the stellar photosphere of a late-type star while the mm emission traces 
a remnant outer disk, but large inner clearings may lead to a lack of mid-IR excess. Better 
coverage of the wavelength range between 10 /im and 1 mm is needed to test this hypothesis. 

For sources without near-IR detections, the "MM" sources in Table [TJ the arguments 
presented above do not apply. Although the emission appears to trace circumstellar dust, 
the fact that no near-IR counterparts are observed suggests high columns of obscuring ma- 
terial. The MM sources all lack mid-IR counterparts as well. These objects appear to be 
so embedded that they are still highly obscured even at 11.7 /im. All of the MM sources 
reside in either the Orion BN/KL or OMC1-S region, bo th of which are known to contain 



young, embedded s ources, HII regions, and outflows (e.g.. IZiurys et al.lll990l ; lBachillerlll996 
Zapata et all 120041 ). 



All of the sources in OMC1-S (LMLA162, MM8, MM13, MM21 and MM22 ) have been 
detected at 1.3 mm wavelength in previous observations (IZapata et al.ll2005l . 120071 ). Measured 
fluxes are similar to those listed in Table [U but somewhat lower in most cases, presumably 
because the poorer uv coverage did not allow large negative sidelobe contributions from 
extended emission to be fully removed . All of these objects a ppear to drive molecular 
outflows traced by CO or SiO emission (IZapata et al.ll2005l . 120061 ). 



We classify the MM sources as candidate Class or Class I protostars. As discussed 
above, it appears that the 1.3 mm emission from at least some of these sources may trace 
dust hotter than 20 K. Such warm dust is expected in the circumstellar environments of 
massive protostars, suggesting that some of the MM sources trace high-mass protostars. 

Several sources detected in previous surveys were not detected here. HC178, HC192, 
HC282, MM3, MM4, MM10, MM15, MM16, M M19, and MM20 should have been detected 
if their 3 mm fluxes (jEisner fc Carpenter! 120061 ) traced dust emission; however, they would 
not have been detected if the objects exhibited flat spectra (e.g., from free-free emission). 
MM7, MM17, and MM18 should have been detected even if they showed flat spectra. We 
detect a 1.3 mm continuum source (MM21) near to HC 178, but find it to be offset by 
more than the relative positional uncertainties, suggesting that the previous association of 
HC 178 with a 3 mm source was mistaken. The other 3 mm objects trace either non-dust, 
potentially time- variable emission, or are spurious, caused by confusion with the BN/KL 
and OMC1-S regions in which they reside. Because our 1 mm observations have far superior 
uv coverage than p revious observat i ons, t hey are less prone to such spurious detections. One 
source detected by I Williams et al.l (120051 ) at 880 /im (171-334) is not detected at 1.3 mm; if 
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the emission comes from small-grained dust, then the expected 1.3 mm flux is comparable 
to our 3a noise level, and hence a non-detection is unsurprising. 



4.2. Frequency of Massive Disks 



We detected 39 sources in our 1 mm mosaics (excluding the BN object). 32 of these 
correspond to (presumed) low-mass near-IR cluster members, and 6 (the "MM" sources) are 
detected only at > 1 mm wavelengths. The remaining detection, LMLA1 62, while not liste d 
as a near-IR source in iHillenbrand fc Carpenterl (120001 ). is seen at 3.6 /im (ILada et al.ll2004j ); 
examination of an archival 2MASS image shows a weak 2 /mi source at this position as 
well. Since the mm-only detections are probably embedded, possibly spherical, protostellar 
objects ( §4.11) . we exclude these from our discussion of disk statistics. Of the remaining 
33 detections, 100% of the 1 mm emission can be attributed to hot gas (free-free) for 6 
sources. Thus, we are left with 27 sources whose 1 mm emission (probably) traces dust in 
protoplanetary disks. 

Since the noise varies across our images, these 27 sources are all detected above slightly 
different thresholds. To examine the frequency of disks more massive than some value, we 
make sensitivity cuts at various levels, examining only the statistics of sources detected above 
chosen noise levels. 

We consider first the 115 cluster members surveyed to a la noise level of 2.7 mJy or 
less. Sources detected above 3a have a circumstellar (dust+gas) mass of > 0.01 M . Nine 
sources (8%) show dust emission of > 8.1 mJy (i.e., 3a detections at this noise level). If we 
use a higher noise cutoff of 5 mJy, then we find that 7 out of 193 stars exhibit dust emission 
above the 3a level of 15 mJy. So, ~ 4% of stars have disks more massive than 0.02 M . If 
we extend the sensitivity cutoff further, to 10 mJy, then 3/254, or ~ 1% of stars are seen 
with disks more massive than 0.04 M . All of the 3a mass levels considered here fall within 



the range of estimates for the MMSN (~ 0.01-0.1 M ; IWeidenschillind Il977l ) . 



The percentage of high-m ass disks derived here can be compared to that determined 
by lEisner fc Carpenterl (120061 ). The observations presented here are substantially more sen- 
sitive than previous observations, and we probe the frequency of disks down to lower mass 
levels; we can there f ore o nly compare statistics for the most massive disks in our sample. 
Eisner fc Carpenterl (120061 ) found that < 2% of cluster members in the ONC have disks more 
massive than 0.1 M . They assumed a distance of 480 pc; their mass limit is actually only 
0.07 M for the distance of 400 pc assumed here. Here we find that < 1% of stars surveyed 
are surrounded by such massive disks, consistent with the estimate from lEisner fc Carpenter 
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( 12006( 1 . 

We emphasize that that results presented above (and in the following sections) depend 
on the conversion of 1.3 mm flux into mass. As discussed in §3.41 there may be some spread 
in the dust properties of our sample that could lead to variations in the derived circumstellar 
masses. For the roughly co-eval, predominantly low-mass cluster population in the ONC, we 
argued that this is a relatively small uncertainty. 



4.3. A Typical Disk in the ONC 



We computed the average flux for the ensemble of non-detected sources in §3.61 The 
average flux indicates that a "typical" non-detected source in the ONC likely possesses a 
disk with a mass of 0.0015 ± 0.0003 M . If we include detected objects (whose dust fluxes 
are listed in Table [2]) in the ensemble, we find that the average disk mass for near-IR cluster 
members in the region is ~ 0.0027 ± 0.0002 M . 

This is comparable to th e average mass determined for 23 proplyds in the ONC at 
880 fim (jWilliams et al.ll2005l ). but substantially lower th an the average mass deter mined 
for > 300 stars at 3 mm wavelengths, 0.005 ± 0.001 M (lEisner fc Carpenter! 120061 ) . The 
discrepancy can be explained in large part by contamination from free-free emission. Because 
the inferred dust mass is proportional to A 3 ^ (Equation [3]), this contamination has a much 
greater affect on the masses inferred from the 3 mm data than on our estimates based on 1 
mm data. The dust mass attributed to free-free emission is (2.3) 3 ~ 12 times larger at 3 mm 
than at 1.3 mm. If free-free emission is present at the < 0.2 mJy level ( §3.31) it would add 



< 0.003 M to the average mass inferred a t 3 mm. Furthermore, if we recompute the mean 
disk mass from lEisner fc Carpenter! (120061 ) using a distance of 400 pc, the average mass is 
decreased by 30%. With the distance correction an d the subtraction of poten tial free- free 
contamination, the recomputed average mass from lEisner fc Carpenter! (120061 ) is > 0.001 



M , in agreement with the estimated average mass inferred from our 1 mm observations. 



4.4. Comparison of Disk and Exoplanet Frequencies 

Less than 10% of stars in the ONC possess disks comparable to the MMSN ( §4.2p . 
Moreover, the average mass measured for the ensemble of (individually) non-detected sources 
is ten times smaller than even the low end of estimates for the MMSN ( §4.3p . indicating that 
the majority of stars do not possess enough mass to form Jupiter-mass planets. 



These statistics can be compared with the frequency of Jupiter-mass planets found 
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around nearby main- sequence stars. 6% of stars surveyed have a Jupiter-mass (or larger) 
planet within 5 AU, while an extrapolation based on current results sug gests up to 10% of 
stars could have a Jupiter-mass planet within 20 AU (IMarcy et all 120051 ). The frequency of 
massive planets is comparable to the frequency of disks in the ONC with (low-end) minimum 
minimum mass solar nebulae. It appears that the MMSN, applied to disk mass measurements 
like those presented here, is a reasonable criteria for forming massive, Jupiter-like planets in 
typical star forming regions like Orion. 



4.5. Disk Evolution 



The frequency of massive disks in the ONC (aged ~ 1 Myr) can be compared with 
surveys of rich clusters of different ages, NGC 2024 (aged ~ 0.3 Myr) and IC 348 (aged 
~ 2 Myr), to constrain the evolution of disks in clustered star forming e nvironments. While 
this c omparison has been made previously using 3 mm observations (e.g.. Eisner fc Carpenter 
20061 ). our 1 mm measurements are less contaminated by free-free emission and yield different 
results ( §4.31) . Although the surveys of NGC 2024 and IC 348 were at 3 mm, the lack of 
O stars in those regions should produce less ionized gas, and hence less contamination by 
free-free emission than in the ONC. 

The average disk masses for "typical" low-mass stars in the three regions is plotted as a 
function of cluster age in Figure [9j For the ONC, we infer a mean disk mass of 0.0027±0.000 2 
In NGC 2024, the mean disk mass is 0.005 ± 0.00 1 M faisner fc Carpenter! I2OO3I ) . 
compared to 0.002 ± 0.001 M in IC 348 dCarpenteri[2002l ). If the differences between NGC 
2024, the ONC, and IC 348 are due to temporal evolution, these observations suggest that 
massive disks/envelopes dissipate on timescales < 1 Myr, and that the average disk mass 
decreases by a factor of 1.9 ± 0.4 between ~ 0.3 and 1 Myr. 



4.6. Dependence of Disk Properties on Environment 



It has been suggested that circumstellar disks in clustered environments may be trun- 
cated due to close encounters with massi ve stars resulting in eit her tidal stripping or photo- 



evaporation of outer disk material (e.g.. IScally &: Clarkd l200ll ). Indeed, photoevaporative 



mass loss has been obse rved from many proply ds, suggesting mass loss rates as high as 



10~ 7 -10- 6 M yr 
of disks over the 



e.g. 



Henney fc Q'Delllll999l ). which would severely deplete the masses 



1 Myr lifetime of the cluster. More detailed models have shown that 
the mass loss rate should be substantially lower for disks with smaller outer radii, since disk 
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material at smaller radii is more tightly gravitationally bound (jClarkdl2007l ). A prediction of 
these models is that larger disks will also be the most massive, since they have to withstand 
higher photoevaporative mass loss rates. 

The p roplyds detected i n our observations are in the top ~ 1/3 of the size distribution 
inferred by lVicente fc Alved (120051 ) . However, the emission from most of these is dominated 
by free-free emission, and even for objects where some component of the flux is due to dust, 
inferred masses are < 0.01 M . Furthermore, there are many other proplyds whose diameters 
are in the top 30% that are not detected in our observations. We also see no obvious trend 
of increasing flux with increasing angular size in our data (Table [2]). Thus, we find little 
evidence that the most extended disks are the most massive. 

Environmental effects on massive disks can also be investigated through the dependence 
of disk properties on cluster radius. We consider the positions of the disks detected in 
our observations (i.e., detected sources corresponding to known near-IR cluster member 
positions) relative to the cluster center, which we define to lie roughly in the middle of the 
four bright Trapezium stars at (a, £)j2ooo — (5 h 35 m 16.34 s , — 5°23'15"6). Figure [TD] shows that 
more massive disks tend to be found further away from the Trapezium stars. If we consider 
only those 194 cluster members where we could have detected disks more massive than 0.02 
M Q above the 3a noise level, we find disks around 1/84 stars (~ 1% ) within 30" and 6/110 
stars (~ 5%) at radii larger than 30". Fisher's exact test indicates 86% probability (1.5a) 
that the small and large cluster radii sources have different frequencies of massive disks. 

Finally, comparison of the ONC with the lower stellar density Taurus region provides 
another test of whether the massive O stars and high stellar density in the Trapezium region 
lead to different disk properties than in more "benign" environments. As discussed in §4.21 
we detected disks more massive than 0.01 M around 9/115 (~ 8%) low-mass ONC cluster 
members. For comparison , 34/153 (~ 22%) of Taurus stars possess such massive disks 
( [Andrews fc Williams! 120051 ) . Fisher's Exact Test yields > 99% probability (3cr) that the 
frequencies of 0.01 M disks in Taurus and the ONC are different. For a slightly higher mass 
cutoff of 0.04 M , such massive disks are found around < 1% of stars in the ONC compared 
to ~ 5% for Taurus. These percentages indicate > 99% probability that the underlying 
distribution of 0.04 M disks in Taurus and the ONC differ. The fraction of approximately 
MMSN-massed disks in Orion is substantially smaller than in Taurus, arguing that the rich 
cluster environment may play a role in limiting the number of massive disks. 



This conclusion differs from that of lEisner fc Carpenter! (12006! ). where the statistics of 
disks more massive than 0.1 M were found to be statistically indistinguis hable in Taurus and 
the ON C. Using the revised distance of 400 pc changes the conclusion from lEisner fc Carpenter 
( 120061 ). because if statistics of (distance-corrected) 0.07 M disks are compared, they are 
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found to be substantially more common in Taurus (Fisher's Exact Test indicates only ~ 1.3% 
probability that the two distributions are the same). Furthermore, the 1 mm observations 
presented here are sensitive to much more of the disk mass distribution, allowing a more 
robust comparison between Taurus and the ONC. 



4.7. Correlation of Circumstellar and Stellar Masses 



Several spectroscopic surveys have provided accurate masses for a subset of the stellar 
popu lation encompassed by our mosaics ([Hillenbrand!! 1 99 71 ; iLuhman et al.ll2000l : ISlesnick et al. 
2004 ). By examining these surveys, after registration of the positions of detected sources to 
the 2MASS grid, we find ~ 130 objects with spectroscopically-determined masses within the 
unit gain contour of our mosaic. While the stellar masses of the remaining cluster members 
contained in our mosaic have been estimated statistically by de-reddening sta rs so that they 
fall on the expected isochrone for the ONC ([Hillenbrand fc Carpenter! 120001 ) the masses of 
individual stars determined in this way have large uncertainties and we do not use them 
here. 

To examine how disk mass depends on stellar mass, we divide the ~ 130 sources with 
spectroscopically-determined masses into three mass bins containing roughly equal numbers 
of stars. The first bin contains stars less massive than 0.3 M , the second bin includes stars 
with masses between 0.3 and 1.0 M , and the third bin contains stars with masses between 
1.0 and 10.0 M . We then make a further cut by excluding all objects for which the noise 
in our 1 mm mosaic at the source position is greater than some cutoff value. As in §4.2[ 
we consider a noise cutoff of 2.7 mJy, which provides a corresponding 3a circumstellar mass 
threshold of 0.01 M . Ideally, we would also bin this sample by clust er radius to contro l 



Hillenbrand! Il997h . 



for potential mass segregation in the inner regions of the ONC (e.g. 
Unfortunately we lack a sufficiently large sample to do this here. 

If we use the raw image fluxes (Tabled]), we find more sources detected in the highest 
stellar mass bin. For 65 stars with spectroscopically-determined masses, surveyed with a 
noise level of 2.7 mJy or lower, we detect 1/21 (5%) stars with M* < 0.3 M , 2/24 (8%) 
stars with 0.3 M < M* < 1 M , and 6/20 (30%) stars with M* > 1 M . However, 
the higher percentage of detected sources in the highest stellar-mass bin is due entirely to 
contamination by free-free emission: if we use dust fluxes from Table [2] then 0/20 stars in 
the highest stellar mass bin are detected. We infer, therefore, that higher mass stars are 
more likely to exhibit free-free emission. 



In contrast, it seems that more massive stars may be less likely to possess massive 
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circumstellar disks. Using the dust fluxes from Table [2J we find that out of the 65 stars 
discussed above, we detect dust emission above the 3a level toward 1/21 (5%) stars with 
M* < 0.3 M , 2/24 (8%) stars with 0.3 M < M t < 1 M , and 0/20 (0%) stars with 
M* > 1 M . This suggests a lower frequency of disks around stars more massive than 1 
M , but these small number statistics do not allow a definite conclusion. A plot of inferred 
circumstellar disk masses versus stellar masses (where available) supports the hypothesis 
that the most massive disks are found around the lowest mass stars. This trend is not, 
however, statistically signi ficant for the small number of objects in Figure [TTJ A similar 
picture was seen in Taurus ([Andrews fc Williams! 120051 ): no correlation between stellar mass 
and disk mass could be established, altho ugh the most mass ive disks were found around 
stars less massive than 1 M . We note that iNatta et al.l (120001 ) claimed to see a correlation 
between disk and stellar masses around early-type stars; however, the dispersion is large and 
the significance of the trend is marginal. 



5. Conclusions 

We imaged a 2' x 2' region of the Orion Nebula cluster in 1.3 mm wavelength continuum 
emission with CARMA and the SMA. Out of > 250 known near-IR cluster members, we 
detected 1.3 mm emission above the 3a noise level toward 33. In addition, we detected 1 mm 
emission above the 5a noise level from six sources not associated with shorter-wavelength 
counterparts. Several of these detected objects are spatially resolved with our observations, 
indicating sizes of ~ 250-450 AU. 

Modeling of long-wavelength fluxes for our targets allowed separation of dust and free- 
free emission components in the measured fluxes. We showed that for the majority of de- 
tected sources, the 1 mm emission appears to trace warm, optically thin, circumstellar dust. 
However, for many of the proplyds, which are located close to the Trapezium stars, the 
millimeter wavelength emission is dominated by thermal free- free emission from hot, ionized 
gas. 

Dust masses inferred for detected sources range from 0.01 to 0.5 M . For the ~ 225 
known near-IR cluster members not detected in our 1 mm observations, images toward the 
positions of near-IR sources were stacked to constrain the mean flux, and circumstellar mass, 
of the ensemble. The average flux is detected at the > 4a confidence level, and implies an 
average disk mass of ~ 0.001 M , approximately an order of magnitude smaller than the 
minimum mass solar nebula. Even when detected sources are included, the average mass is 
< 0.003 M . While the derived masses are uncertain by a factor of 3 or so (mostly due to 
uncertainties in the dust opacity), the range of possible average disk masses is still smaller 



-23- 



than even the low end of estimates for the MMSN. A "typical" star in the ONC does not 
appear to possess sufficient mass in small dust grains to form Jupiter-mass (or larger) planets. 
Evidently, giant planet formation is either advanced (having thus depleted the small dust 
grains in the disk) or impossible around most stars in the ONC. 

We compared the average disk mass inferred for the ONC with similarly determined 
average masses in older and younger clusters. We find evidence for evolution of the dust 
(most likely depletion or agglomeration) on ~ 1 Myr timescales. Between ~ 0.3 Myr and 
~ 1 Myr, the average disk mass decreases by a factor of 1.9 ± 0.4. 

The percentage of stars in Orion surrounded by disks more massive than the minimum 
mass solar nebula is < 10%. This is significantly lower than in Taurus, indicating that en- 
vironment has an impact on the disk mass distribution. Our data suggest (with marginal 
statistical significance) that the most massive disks may be located further from the Trapez- 
ium stars, supporting the hypothesis that photoevaporation may be truncating disks near 
to the cluster center. 

Finally, our observations show no clear correlation between stellar mass and disk mass, 
but suggest that massive disks may be more likely to be found around lower mass stars. 
The percentage of detected disks is lower for stars more massive than 1 M , and the most 
massive disks detected are associated with the relatively low stellar mass stars in the sample. 
However, larger numbers of stellar and disk mass measurements in the ONC are needed to 
build up better statistics and further constrain the relationship between stellar and disk 
properties. 
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A. Free-free emission spectrum for an ionized wind 

We compute the spectrum expected for free-free emission from ionized gas whose density 
depends on stello centric radius as 

n gas = n (R/Rv)- a . (Al) 

For an ionized wind with constant M, such as we expect for proplyds, a = 2. However, 
massive stars may ionize their circumstellar environments directly, in which case the density 
distribution may differ from this R~ 2 power-law. 



For a power-law density profile (Equation IA1I) . the optical depth can be approximated 



as (e.g., 



Altenhoff et al.lll960f ). 



^ r 0.16ng dz 

Here, z is the line of sight through the ionized gas in pc, x is the impact parameter in pc, no 
is the normalization of the gas density in cm -3 (assumed to represent the electron and ion 
densities), T e is the electron temperature, and v is the frequency in GHz. This expression 
can be integrated directly as long as a > 1.5, with the result 

^ O.OSngi* _ (2a _ 1} ^(-1)^(1/2) 
U ' S ~ z/ 2 - 1 ^ 1 - 35 sin(7r/2)(a-l)!r(l/(2-a + l))' 1 ; 

The last fraction has order unity, and we can thus approximate r as 

T ^ ff ~ ,.2.1 T 1.35 X ■ ^ A4 J 

e 

We can now invert Equation IA4I to determine the maximum impact parameter for which 

'0.08n 2 R^ 



- ( ) ^ (A5) 



For impact parameters larger than x T ~i, the gas is optically thin, and at smaller impact 
parameters the gas is optically thick. The total spectrum of free-free emission from the 
source can be approximated by the blackbody flux times the solid angle of the optically 
thick region, as long as the optically thick region is finite. When the entire wind becomes 
optically thin, the spectrum flattens. 

We can thus parameterize the free-free emission from a wind as 

-^.turn^/t'turn) -0 ' 1 if V > U tuln 

"■'i- " : S P / / (A6) 
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see 



Wright fc Barlowlll975l . for a somewhat different derivation of this result). For a spher- 



ical wind with constant M, a = 2 and the optically thick material emits with F v oc z/ 6 . 
Steeper exponents can result if the ionized gas density drops off more steeply than R~ 2 , 
for e xample as might occu r in a centrally illuminated wind from a massive star like BN 
Plambeck et al.l[l995h . Alternatively, non-spherical wind geometries can lead to shal- 



le.g. 



l ower radial density pro files, and hence shallower spectral slopes of the free-free emission 
rtWright fc Barlow! Il975h . " ? or simplicity, since we generally have a limited number of data 
points with which to constrain the free-free emission spectrum, we will assume a = 2 in the 
present analysis. 
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Table 1. Long- wavelength fluxes for sources detected in A1.3 mm continuum 
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Table 1 — Continued 



ID 


a 

(J2000) 


5 

(J2000) 


(mjy) 


5*1 . 3mm 

(mjy) 


53 mm 

(mjy) 


6mm 

(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


MM13 
MM21 
MM22 
MM23 
MM24 


5 35 13.75 
5 35 13.57 
5 35 13.66 
5 35 14.00 
5 35 14.62 


-5 24 7.74 
-5 23 59.04 
-5 23 54.94 
-5 22 45.04 
-5 22 28.94 




317.5 ±25.1 

153.6 ± 13.2 
123.0 ± 12.7 

61.2 ±7.4 
167.4 ±27.0 


36.9 ±8.3 
15.6 ±4.7 
17.1 ±4.2 
9.6 ±3.1 
< 32.7 




0.9 ±0.1 
0.4 ±0.1 









Note. — All quoted uncertainties are la; uncertainties (of 10%) in the overall flux scale are not includ ed in this table. The MM 
source s are detected only at wavelengths > 1 mm. HC objects are known near-IR cluster members from Hillenbr and fe Carpenterl 
d2000h that are also detected at m m wavelengths, and the sources with numerica l labels are a subset of th i s sam ple that are also 
detect ed optically as pro plyds (e.g.. [O^Dell et al]|l993l ). LMLA162 is absent from [Hillenbrand fe Carpenterl JgOOO), but seen at 3.6 
bvlLada et al.l [2 004): we also see a faint 2 fim point sou rce at this location in a 2M ASS image. MM8 and MM13 were detected 
bv lEisner fc CarpenteTil2006h . and MM21 was detected by lEisner fc Carpente 2 J2006h but associated with HC 178 (we find that 
MM21 and HC 178 are separated by an angle larger than our rel ative position a l unc ertainties). LMLA162, MM8, MM13, MM21, 

and MM22 were also detected in previous 1.3 mm observations bv lZapata et al] ll2005h . MM23 and MM24 are newly detected here. i — > 

References. — References for the fluxes arc as follows: 880 um [IWilliams et al.l |2005|); 1.3 mm (t his wor k) ; 3 mm 
( Eisne r fc Carpente r 2006): 3 .6 mm llMundv et aljl995h : 1.3 cm dZapata et alj2004) ; 2 and 6 cm et al ■Il993bl) ; 20 cm jFelli et alj 

Il993al) . Quoted upper limits for the 3 mm fluxes are 3<r. 
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Table 2. Derived quantities for detected sources 



ID 


S^.dust (mJy) 


A'-fcircumstellar (Mq) 


R dia k (AU) 


s v,t>\ ( mJ y) 


147-323 


7.5 ± 2.2 


0.009 ± 0.003 


88 


83 


155-338 


4.0 ± 4.0 


0.005 ± 0.005 


102 


112 


158-323 


1.0 ± 1.0 


0.001 ± 0.001 


105 


119 


158-327 


5.0 ± 4.0 


0.006 ± 0.005 


122 


161 


159-350 


33.0 ± 7.0 


0.042 ± 0.009 


152 


250 


163-317 


6.0 ± 2.0 


0.008 ± 0.003 


93 


93 


167-317 


2.0 ± 2.0 


0.003 ± 0.003 


122 


161 


168-326NS 


7.0 ± 4.0 


0.009 ± 0.005 


<100 


< 108 


170-337 


10.0 ± 2.5 


0.013 ± 0.003 


126 


171 


171-340 


13.0 ± 2.3 


0.016 ± 0.003 


80 


69 


177-341 


4.0 ± 4.0 


0.005 ± 0.005 


177 


339 


HC180 


10.7 ±3.0 


0.013 ± 0.004 


<100 


< 108 


HC189 


99.5 ± 8.4 


0.125 ± 0.011 


<100 


< 108 


HC246 


17.8 ± 2.4 


0.022 ± 0.003 


<100 


< 108 


HC254 


17.7 ±3.8 


0.022 ± 0.005 


<100 


< 108 


HC295 


8.0 ± 3.5 


0.010 ± 0.004 


187 


381 


HC336 


0.0 ± 0.0 


0.000 ± 0.000 


<100 


< 108 


HC350 


5.0 ± 3.5 


0.006 ± 0.004 


<100 


< 108 


HC351 


8.7 ± 2.8 


0.011 ± 0.003 


<100 


< 108 


HC361 


19.2 ± 3.2 


0.024 ± 0.004 


<100 


< 108 


HC383 


7.0 ± 2.3 


0.009 ± 0.003 


202 


111 


HC401 


7.3 ± 2.1 


0.009 ± 0.003 


<100 


< 108 


HC412 


7.7 ± 2.3 


0.010 ± 0.003 


161 


283 


HC414 


3.0 ± 3.0 


0.004 ± 0.004 


<100 


< 108 


HC418 


7.6 ±2.5 


0.010 ± 0.003 


<100 


< 108 


HC436 


1.0 ± 1.0 


0.001 ± 0.001 


<100 


< 108 


HC438 


67.8 ± 14.2 


0.085 ± 0.018 


200 


433 


HC440 


13.0 ±2.8 


0.016 ± 0.004 


<100 


< 108 


HC495 


22.1 ±6.6 


0.028 ± 0.008 


<100 


< 108 


HC498 


15.3 ±4.6 


0.019 ± 0.006 


<100 


< 108 


HC514 


23.7 ± 7.4 


0.030 ± 0.009 


166 


299 


HC771 


16.1 ±4.9 


0.020 ± 0.006 


176 


337 


LMLA162 


103.3 ± 7.2 


0.130 ± 0.009 


<100 


< 108 


MM8 


407.2 ±27.5 


0.512 ± 0.035 


<100 


< 108 


MM13 


317.5 ±25.1 


0.399 ± 0.032 


<100 


< 108 


MM21 


153.6 ± 13.2 


0.193 ± 0.017 


<100 


< 108 


MM22 


123.0 ± 12.7 


0.155 ± 0.016 


<100 


< 108 


MM23 


61.2 ±7.4 


0.077 ± 0.009 


<100 


< 108 


MM24 


167.3 ±27.0 


0.211 ± 0.034 


217 


514 



Note. — S^dust is the component of the observed 1.3 mm emission due to cool dust, 
determined from a fit to long-wavclcngth fluxes of a model including thermal free-free 
emission as well as dust emission (see Figure [6j. Quoted uncertainties are la, and 
do not include systematic uncertainties associated with the overall flux calibration or 
conversion from flux to mass. 
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Fig. 1. — Pointing positions for the CARMA mosaic (triangle symbols) and the SMA mosaic 
(square symbols) plotted over the positions of if -band sources in the ONC (op en circles). 
These if-band source positions are drawn from Hillenbrand Carpenter! (120001 ). and have 
been registered to the 2MASS astrometric grid. The unit gain contour of the CARMA 
mosaic (dotted curve), the SMA mosaic (dashed curve), and the mosaic produced from the 
combined SMA+CARMA dataset (solid curve) are also indicated. 
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Fig. 2. — The Orion Nebula cluster, imaged in A1.3 mm continuum with CARMA. Only 
data observed on long baselines (r uv > 70 kA) were used to create this image. The unit 
gain region of the mosaic is indicated by a solid curve. The synthesized beam has FWHM 
of 0''61 x 0"52 at a PA of 70°. The RMS noise varies within the unit gain contour from 2.3 
to 23 mJy (la), with a mean RMS of 4.6 mJy Partially resolved extended emission in the 
BN/KL and OMC1-S regions increases the noise level in the upper and lower right corners 
of the image. The map has been divided by the theoretical RMS, to visually down-weight 
the noisier edges. 
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Fig. 3. — The Orion Nebula cluster, imaged in A1.3 mm continuum with the SMA. Only 
data observed on long baselines (r uv > 70 kA) were used to create the image. The unit 
gain region of the mosaic is indicated by a solid curve. The synthesized beam has a roughly 
circularly symmetric core with a FWHM of 0"98. The RMS within the unit gain contour 
varies from 0.8 to 28 mJy with a mean value of 2.7 mJy The BN/KL and OMC1-S regions 
produce the strong emission visible in the upper and lower right corners of the map. The 
map has been divided by the theoretical RMS, to visually down- weight the noisier edges. 
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Fig. 4. — The Orion Nebula cluster, imaged in A 1.3 mm continuum with CARMA and 
the SMA (greyscale). CARMA and SMA data were combined in the uv plane, and only 
data observed on long baselines (r uv > 70 kA) were used to create this image. The angular 
resolution is ~ 0'.'7 x 0''6. The unit gain region of the mosaic encompasses a 2' x 2' area, 
as indicated by the solid contour, and the RMS residuals (lcr) within this region vary from 
1.8 mJy in regions devoid of bright emission to > 10 mJy in the crowded regions toward the 
right of the map. Millimeter detections above the 3cr level coincident with near-IR cluster 
members and proplyds, and sources without counterparts detected above the 5a level ( §3. lj) . 
are indicated by gray circles. The map has been divided by the theoretical RMS, to visually 
down-weight the noisier edges. 
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Fig. 5. — Contour images of sources detected in 1 mm continuum emission. Contour incre- 
ments are la, beginning at ±2<r, where a is determined locally for each object (see §3.ip . 
Solid contours represent positive emission, and dotted contours trace negative features. The 
BN object is excluded from these plots, as it is discussed in more detail in a later paper. For 
sources detected at infrared wavelengths, individual images are centered on the previously 
measured near-IR coordinates. For the MM sources, detected only at > 1 mm wavelengths, 
images are centered on the peak fluxes. 
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Fig. 5. - 



continued. 




Fig. 6. — Millimeter and radio fluxes for our sample (points), along with best-fit models 
including free-free and thermal dust emission. Models including free-free and dust emission 
are indicated by solid lines, and dotted and dashed lines show dust-only models with (5 — 
and /3 — 1, respectively. The free-free flux density is proportional to z/~ ai for gas that 
is optically thin for all radii, and to u 0M for a partially optically thick gaseous wind; the 
emission is thus parameterized by the flux at a single wavelength and a turnover frequency. 
The emission from cool dust is proportional to z/ 2+/3 . 
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Fig. 6. - 



continued. 
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Fig. 7. — Average image, obtained by stacking the 1 mm continuum emission observed 
toward each of 226 low-mass near-IR sources not detected individually above the 3a level. 
Contour levels begin at ±2<r and the contour interval is la (negative contours are shown as 
dotted lines). Emission is detected for the ensemble at a significance of > 4a, and exhibits 
a compact (and beam-like) morphology approximately centered on the origin. The degree 
to which the emission is smeared out is consistent with the ~ 0''4 positional errors in the 
near-IR source positions. 
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Fig. 8. — Contour image of the 1.3 mm wavelength continuum emission toward the proplyd 
177-341, overlaid on an HST map (IBally et al.ll2000l ) of the same position. Positive contours 
are white, solid curves and negative contours are dotted black curves. Contours are ±2, 4 
and 6a. The FWHM of the synthesized beam is shown as a filled back symbol. 
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Fig. 9. — Average disk mass as a function of age for the NGC 2024, O NC, a nd IC 348 c l usters . 



The disk masses are taken from this work, Eisner &: Carpenter 



(20031). an dlCarpente r (2002) 



(1997 


), 


Luhman et al. 


(1998 


), and 


Luhman 


(1999) 



IC 348 were measured at 3 mm, where potential contributions from free-free emission would 
be stronger than for the average mass measured here for the ONC at 1 mm. We argue in 
751 however, that free- free contamination is unlikely in NGC 2024 and IC 348. 
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Fig. 10. — Circumstellar mass as a function of radial distance from the center of the Trapez- 
ium stars for detected objects. 
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Fig. 11. — Circumstellar mass as a function of stellar mass for the subset of detected 
objects where spectroscopically determi ned stellar masses are available ([Hillenbrand! 119971 ; 
Luhman et aDbood ; Isiesnick et aPhood ). 



